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We study solitary-wave and kink-wave solutions of a modified Boussinesq equation through a 
multiple-time reductive perturbation method. We use appropriated modified Korteweg-de Vries 
hierarchies to eliminate secular producing terms in each order of the perturbative scheme. We 
show that the multiple-time variables needed to obtain a regular perturbative series are completely 
determined by the associated linear theory in the case of a solitary-wave solution, but requires 
the knowledge of each order of the perturbative series in the case of a kink-wave solution. These 
appropriate multiple-time variables allow us to show that the solitary-wave as well as the kink- 
wave solutions of the modified Botussinesq equation are actually respectively a solitary-wave and a 
kink-wave satisfying all the equations of suitable modified Korteweg-de Vries hierarchies. 



I. INTRODUCTION 



In this paper we study two model equations of the modified Boussinesq type 

u u - u xx + u xxxx + 2a{u p ) xx = , (1) 

with u{x,i) a one-dimensional field, a = ±1, p = 3 and subscripts denoting partial differentiation. These equations 
approximately describe the propagation of waves in certain nonlinear dispersive systems. For example they appear 
when we study the continuous limit of a Fermi-Pasta-Ulam dynamical system with cubic nonlinearity (Dodd et al. 
1982). They can appear also governing the evolution of long internal waves of moderate amplitude, (Ablowitz et al. 
1973) or describing the dynamic of a stretched string (Ablowitz et al. 1973; Mott 1973 ). 

They are considered as intermediate long-wave equations since they represent an intermediate dynamic, in com- 
plexity and completeness, situated between the complete dynamic of the full initial equations describing any wave 
number and any amplitude, and some strong long-waves and small- amplitude limits. Their appropriated long-wave 
and small-amplitude limit with a further restriction to unidirectional propagation yields the modified Korteweg-de 
Vries equation (mKdV). 

Asymptotic methods are very often employed and very useful to study problems of this type. Hence, to study the 
evolution of long-waves of equations (0), we will consider a perturbative scheme based on the reductive perturbation 
method of Taniuti (Taniuti 1974). We will include the classical slow space-variable £, but we will modify it by 
introducing an infinite number of slow time-variables: T3, T5, T7, etc. As we are going to see, by using these slow 
variables together with all the equations of an appropriated modified Korteweg-de Vries hierarchy (Chern & Peng 
1979 ) , we can obtain the solitary-wave solution of the modified Boussinesq equation (Q) with a = 1 (for shorteness 
mBI) 

u(x, t) = k sech jfc (x - \/\ - k 2 ^j tj , (2) 

as well as the kink- wave solution of the modified Boussinesq equation (JlJ) with a — — 1 (mBII) 

u(x, t) = k tanh |fc (x - \J \ + 2k 2 ^j tX . (3) 

The solutions @) and (|^), written in the laboratory coordinates (x,y), are actually built up respectively from the 
solitary-wave solution and the kink-wave solution of the whole set of equations of appropriated mKdV hierarchies 
written in the slow space- variable £ and in each one of the slow time-variables T3, T5, ry, ... 

These results follow from the general long-wave perturbation theory and from the requirement of uniformity for 
large time of the associated perturbative series. This last fact makes the perturbative series truncates for solutions 
of type (|^) or (|3|), rendering thus exact solutions for the modified Boussinesq equations. Furthermore, we will show 
that the elimination of the secular producing terms in the perturbative series linked to the solution (|J), is completely 
accounted for only the linear theory associated to M), which allows us to know "a priori" all the constants which 
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define the slow times-variables r s in function of t. These properly normalized slow time coordinates automatically 
give us a perturbative series which is free of secularities (uniform expansion). 

On the other hand, to solve the same problem for the solution (||), that is, to eliminate the secular producing terms 
linked to solution (|§|), requires the knowledge of each term of the perturbative series. 

This paper is organized as follows. In Sec. II the multiple time formalism is introduced for the modified Boussinesq 
equations (|l|), the first few evolution equations are obtained and the problems associated with them are exhibited. In 
Sec. Ill we discuss how the mKdV hierarchy shows up. In Sec. IV we show how they can be used, in the case a = 1, 
to eliminate the soliton related secularities of the evolution equations for the higher-order terms of the wave fields. 
In Sec.V by returning from the slow variables to the laboratory coordinates we obtain the above mentioned relation 
between the solitary-wave of mBI and the corresponding modified Korteweg-de Vries hierarchy. Sec. VI is consecrated 
to the case of mBII (a = — 1). In Sec. VII we give a general proof that the symmetries of time derivatives lead to the 
modified Korteweg-de Vries hierarchy. Finally in Sec. VIII we summarize and discuss the results obtained. 

II. THE MULTIPLE TIME FORMALISME: MIXED-SECULAR TERMS 

In order to study the far-field dynamics of long-wave solutions of eq. ([!]), we will need to define slow space and time 
variables. A small parameter e giving the long-wave character of the studied solution is introduced via the definition 

k = en, (4) 

where k is the wave number and k a parameter of order one. Accordingly, we define a slow space variable 

e = <x ~ t) , (5) 

as well as an infinity of slow time coordinates 

T 3 = e 3 t, T 5 = e 5 t, T 7 = e 7 t,.... (6) 

Consequently, we have 

dd~drd 7 d 

dt = - e di + e d^ + e d^ + e Wr + --" (8) 

Moreover, we consider a small-amplitude solution of (|l|) and we make the expansion 

u = eu = e(uo + e 2 u 2 + e 4 U4 + ...). (9) 

We suppose that U2n = u 2n (^, 73, T5, ...) for n = 0,1,2,... which corresponds to an extension of the function u in 
the Sandri's sense (Sandri 1965). Substituting eqs (0) (JsJ) and (||) into the Boussinesq equations (0) the resulting 
expression, up to terms of order e 4 , is: 

d 2 d i 2 fd 2 n d 2 \ 4 / „ d 2 n d 2 
ttt. + — r + e 2 — ^ - 2—— + e 4 -2—— + 2— — 



and 



_2 1 h e' 1 - 1 > 

<9£<9t 3 <9£ 4 \9t 2 <9£<9t 5 / \ d£,dr 7 dr 3 dr 5 

d 2 

2a ^H N + e 2 (3u 2 u 2 ) + e\Zu 2 Ui + 3u Q u 2 ) + ...] = , (10) 



At order e , after an integration in £, we get 



du d 3 u 2 du 



which is the mKdV equation. 

At order e 2 , eq.din) yields, using (|ll|) and integrating once in £ 
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du 2 , d 3 u 2 d(v%u 2 ) du ld 5 u /<9u \ 3 n du d 2 u 

2 d 3 u 2 4 du 



" 3aM o^F-9o '<~ojr- (12) 

Equation ( |l2| ) is a linearized inhomogeneous mKdV whose general solution consists of a sum of a general solution to 
the homogeneous equation and a particular solution to the nonhomogeneous equation. As it stands, (|T^), presents two 
problems. First, the inhomogenity (source term) is unknown because the evolution of uq in the time T5 is not known. 
The second problem is related to nonuniformity of the expansion for u. When we considerer a soliton type solution of 
eq.([ll]), case a — 1, the term Q^f- is proportional to which is a solution of the associated homogeneous equation. 
Hence the general solution u 2 of ([[2]) contains a term proportional to T3 which gives rise to a nonuniformity in 
the perturbative series (mixed-secular term). 

For mBII, case a = —1, we have kink-type solutions to eq. (O). Solutions of type (||) exist but are complex 
and we will not consider them. In this case the linear term of eq.(|12[), and actually some nonlinear ones, produces 
secularity as well. 

In the next sections we will deal with these two problems. 

III. THE RISE OF THE MODIFIED KORTEWEG-DE VRIES HIERARCHY 

As we have seen, the field uq satisfies the mKdV equation in the time T3. The evolution of the same field uq in any 
of the higher order times T 2n +i can be obtained in the following way (Kraenkel et al. 1995). 

First, to have a well ordered perturbative scheme we impose that each one of the equations for mo,t 2 „+i be e- 
independent when passing from the slow variables (uq r 2n+ i) to the laboratory coordinates (u,x,t). This step 
selects all possible terms which can appear in wo,r 2 „ +1 - For instance, the evolution of uo in T5 is restricted to the form 

u ,r 5 = au 0M + bulu 0M + cu uo^uo i2 ^ + d(u ,^) 3 + eu^uo,^ + fu u 0AS 

+ guluo,2t; + huQ^UQ^ + iul + j(u 0i2 ^) 3 , (13) 

where a, b, c, d, e, f, g, h, i and j are unknown constants. Then, by imposing the natural (in the multiple time 
formalism) compatibility condition 

( U Q.T 3 ) = (uo,T 2n+1 ) , (14) 

with uq satisfying mKdV in T3, it is possible to determine any uo, T2 „ + i i-e to determine all the contants appearing in 

M 0,r 2 , 1 + i- 

As it will be shown in Sec. VI, the resulting equations are those of the mKdV hierarchy. In particular, for uo. T5 and 
mo,t 7 , using the mKdV in its canonical form (for shorteness mKdVI) with a — 1 

Uo,t 3 + wo,3£ + GuIuqx = 0, (15) 

we obtain 

uo,r 5 = a 5 ruo,5£ + 10uqU ,3£ + 40u uo,^uo,2C + 10(u ,c) 3 + 30uqM ,c) , (16) 

and 

w ,r 7 = -ar( + Uo,7(, + 84m wo,4Uo,45 + 560woU ,£Wo,2C + 14uoWo,5? + 140u w , 24^0,35 
+ 126u ,3 5 Ke) 2 + 182u , ? K2 5 ) 2 + 70(u ) 4 u M + 420(u ) 2 (wo,«) 3 
+ UOuo^uo) 6 ) . (17) 

For the alternative form, mKdVII (a = —1) 
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u ,t 3 + "0,3? ~ 6wqU 0j? = , (18) 

i*o,r 5 and uo,t 7 are given by 

uo,r 5 = Pa(uoM ~ 10moUo,35 - ^OuqUq.^uq^ - 10(u , 4 ) 3 + 30uqU ^ , (19) 

uq,t 7 =/3 7 (- it , 7£ + 84u uo,eUo,45 - 560uoi*o,^u ,25 + IAuqUq^ + 140u i*o,2£Uo,3£ 
+ 126i*o,3cK,«) 2 + 182i*o, 4 (uo,2e) 2 - 70(it ) 4 i*o,3C - 420(u ) 2 (i*o,«) 3 
+ 140u , C K) 6 ) • (20) 

The coefficients as, a?, 05, (3-?, are free parameters which are not determined by the algebraic system originated from 
eq. (Q) . These free parameters are related to different possible normalizations of the slow time variables and we will 
use them to obtain a secular free perturbative solution of ([!]). 

We will see in Sec.V and VI that the genesis of the values of these parameters are rather different for the perturbation 
theory associated with the two cases a = 1 or a = - 1. 

IV. HIGHER ORDER EVOLUTION EQUATIONS FOR MODIFIED BOUSSINESQ I 

Let us consider the perturbation theory asociated with mBI (a = 1 in eq.((lj)) with T5, tj, as defined by (g) and 
T3 redefined by 

r 3 = (21) 

Hence the operator (||) now reads 

d d lo<9 . <9 7 d . „. 

wt = -%-rd7 3 +e d^ +e d7 7 + --- (22) 

At first order, this definition of T3 gives the mKdVI in canonical form (|l^). At order e 2 we obtain the following 
equation: 

1 /■« 

£i(ii 2 ) = 2uo,r 5 - 7 / u o,2r 3 ^ , (23) 

where Cj is the linearized mKdVI operator defined by 

Cj(v)=v T3 +v 36 +6(ulv) 6 . (24) 
Substituting i*o,r 3 and i*o,r B from ( |l5[) and (|l6| ) we have 

Ci{u 2 ) = 2(a 5 - ^j u 0M + ( 20q; 5 - 3^iiQit ,3C + (sOcts - 9^u i*o,4"o,2C 

+ (20a 5 - 3) (uo.t f + (60a 5 - Qji^uo^. (25) 

If we assume the solution of the mKdVI (eq.([l5|)) to be the solitary- wave solution 

u = k sech (k£ — k 3 t 3 + 6), (26) 



where 9 is a phase, it is easy to see that only the term 1*0,55 ^ s a secular-producing term. It can be eliminated if we 
choose as = A. In this case, eq.dga) becomes 
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Ci{u 2 ) = --u 2 u ^ + u Q u Q . ( u ^ - -^(uq.^) 3 - 2 U o w o,«- (27) 



Using for uq the solitary- wave solution (Eq), we see that the right-hand side of eq. (|27|) vanishes, leading to 
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Cx{u 2 ) = 0, (28) 
which is a homogeneous linearized mKdVI equation, for which we will consider the trivial solution 

u 2 = 0. (29) 
At order e 4 , and already assuming that 1*2 = 0, we obtain 

Cx(u 4 ) = 2u ,r 7 + / Uo,r 3 ,r 5 rfC • (30) 
J — 00 

Using equations (15) and (|l6|), with = |, to express respectively u T3 and u T5 we obtain 

1 \ 65 „ 35 



/ ! \ bo 00 



4? 



O 2 « 2 105 2 2 7 , 3 45 4 

- 2 u ,55M - y u o,5(wo,2c) 2~w (u ^) - 75w ,£«0,2£ w ~ y u o,3|M 

-fW(n ) 6 . (31) 

Now the source term proportional to tio,7£ is the only resonant. Then, in the same way we did before, we first use 
the higher mKdVI, eq.(|l7), to express uq jT7 and then we choose the free parameter a? in such a way to eliminate the 
resonant term uqj^. This choice corresponds to 07 = - jj, which brings eq. (|3l|) to the form 

5 

+ jUoA u o,2t;) 2 - 2 U $M u o ~ 5 "o,£«o ■ (32) 



Substituting again the solitary-wave solution (g6|) for uq, we see that the nonhomogeneous term of eq.(|32j) vanishes, 
leading to 

C T {u 4 ) = 0. (33) 

Again, we take the trivial solution 

m 4 = 0. (34) 

This is a general result that will repeat at any higher order. For n > 1 the evolution equation for u 2 n, after using 
the mKdVI hierarchy equations to express uo tT2n+1 and after eliminating the secular producing term coming from the 
solitary-wave solution (|26]), is given by an homogeneous linearized mKdVI equation. 
Consequently, the solution U2n = 0, for n > 1, can be assumed for any higher order. 



V. BACK TO THE LABORATORY COORDINATES. CONNECTION TO THE DISPERSION RELATION 

Actually, to obtain a perturbative scheme free of secular producing terms at any higher order, we assume that uq 
is the solitary-wave solution of all the equations of the mKdVI hierarchy, each one in a different slow-time variable. 
Such a solution may be obtained and is given by 



Uo(£,T3,t 5 ,T7, ...) = k sech 



ft£ — k 3 t 3 + n 5 a 5 T 5 — n 7 a 7 T 7 + .. 



(35) 



First, recall that we have expanded u according to eq.(^J). Thereafter, we have found a particular solution in which 
U2n — for n > 1. Consequently expansion (M) truncates, leading to an exact solution of the form 



u = euo , 



(36) 



with Mo given by eq. (|35|) . Moreover, the slow variables (k, £, T2 n +i) are related to the laboratory ones ( k, x, t ), 
respectively by eqs. (|4J), (0), (|^) and ( pl| ) and we have found later that 



a5 = 8' a7 = "l6' 



(37) 
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Then, in the laboratory coordinates, the exact solution ( |36| ) is written 

u(k, x,t) — k sechA; 



_ ( 1 - -k 2 — -k A — —A; 6 - 
2 8 16 



Now, the series appearing inside the parenthesis can be summed: 



1 _ Ifc 2 - ifc 4 - —A; 6 
2 8 16 



Vl - k 2 



Consequently, we get 



k sech k (x — \fl — k 2 tj . 



which is the solitary-wave solution of mBI. 

Let us now take the linear mBI dispersion relation 



w(Jfc) = k y/l + k 2 . 



Its long-wave expansion (k = en) is given by 



uj(k) = en + l-e 3 K a - ^e 5 K 5 
2 8 



16 



7 7 

e k 



(38) 



(39) 



(40) 



(41) 



(42) 



In passing we notice that the absolute value of the coefficients of the expansion coincide exactly with those found: 
a3 to obtain the mKdVI in canonical form and «5, otj, necessary to eliminate the secular producing term in each 
order of the perturbative scheme. 

With this expansion, the solution of the associated linear mBI reads 



u = exp i 



:{ X -t) 



i«[.r-l\ - -e 3 K?t + g£ 5 K 5 *- jj; 



e 7 K 7 t + 



(43) 



Therefore, if we define from the begining, as given by this expression, the properly normalized slow time coordinates 

(44) 
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1 3 1 5 1 7 



the mKdVI (15) will be obtained at first order and the resulting perturbative theory will be automatically free of 
secular producing terms for the soliton solution. 



VI. HIGHER ORDER EVOLUTION EQUATION FOR THE MODIFIED 

BOUSSINESQ II 

We consider now the case a = — 1 in eq.(Q). We assume the same operators and given by (0) and (p2|), and 
the same expansion for u. We obtain at order e° and e 2 the following equations 

«o,T 3 + w ,3C - 6uo^o,c = 0, (45) 

which is mKdVII, and 
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Cji{u 2 ) = 2u ,T 5 - 4 / uo,2r 3 d£, , (46) 



where Cxx is the linearized mKdVII operator. Using eq. (|45[) to express uq iT3 and eq.(19) to express wo,t 5 in (|46j) we 
obtain 

Cxx(u 2 ) = (2(3 5 - ^ju 0y5S - (20/3 5 - 3^UqM ,3C - (80/3 5 - ^uqUq^uq^ 

- (20/35 - 3) K ? ) 3 + (60/3 5 - 9) u%u QA . (47) 
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In this case the real solution of mKdVII is of kink-wave type. It reads 

u (x,t) = k tanh (k£ + 2k 3 t 3 + 9) . (48) 

When we substitute ( ft8[ ) in (f47j), resonant terms appear. Contrary to the case of mBI, they do not come only from 
the linear term, but also from the nonlinear ones. They are proportional to 



u ,£ = k secli Ik£ + 2k j t 3 + 6j , (49) 

and we have 

Cix{u 2 ) = -K sech 2 (k£ + 2k 3 t 3 + 9\ (l - 12/3 5 ) . (50) 

Hence, the value 

A - T2 < 51 > 

eliminates the secular producing term at order e 2 , and gives Cxziuv) = 0. We then assume the trivial solution u 2 = 0. 
At order e 4 we have 

Ctz(v>4,) = 2m 0:T7 + / u ,T 3 ,T 5 d^ . (52) 



Using dl8[ ) ( |19| ) and ( |20| ) to express respectively Mo,r 3 , M o,r 5 and u ,r 7 we obtain, (using (|4q)) 

£zzM = « sech 2 + 2k 3 t 3 + 6 1 ) (l + 40/3 7 ) . (53) 

We choose 

(h = (54) 
to eliminate the secular producing term at this order. Again, we take the trivial solution 

u 4 = 0. (55) 
Actually, as in the previous case, we can assume U2n = for n > 1 and we have the exact solution 

it = e 2 uo, (56) 
where uo is the kink-wave solution to all equations of the mKdVII hierarchy, which reads 

u = Ktahn (k£ + 2k 3 t 3 + 6k 5 /3 5 t 5 + 20k 7 /3 7 t 7 + . (57) 



In the laboratory coordinates, the exact solution (|57|) is written as 

I. a 1 



u 



k tahn k 



a;- (l + fc 2 - -fc 4 + -fc 6 -...)<]. (58) 



The series appearing inside the parenthesis can be summed: 



1 + k 2 - ifc 4 + ifc 6 - ... = v 7 ! + 2fc 2 , (59) 



and we obtain the kink-wave solution of mBII 



u = k tanh [kx - ky/l + 2k 2 tj . (60) 

To obtain (^0|) we used the definitions 

ra = ~\e\ r 5 = lA, r 7 = -l e *t, ... . (61) 

In this case the cofficients are not those obtained from the long-wave expansion of the linear dispersion relation of 
mBII. 



VII. COMMUTATIVITY OF TIME DERIVATIVES IMPLIES THE MKDV HIERARCHY: A GENERAL 

PROOF 

In this section we give a general proof that the symmetries of time derivatives, with the scale invariance requirement, 
leads to the mKdV hierarchy in T5, Tj,..., for a field which satisfies mKdV in T3. We do this for mKdVI. 
Let M n be a polynomial in terms of the form 

u a °ufull...u a ^ (62) 

with a,; e N. We define the Rank (R) of M n as 

l 

fl(M„)=^(l+j) % . (63) 
3=0 

Let us now consider our principle which requires that each one of the evolutions equations 

u T3 , u T5 , u T7 ,...,u T2n _ 3 n = 2, 3, (64) 

be e- independent when passing from the slow (uo>£> ^n-i) to the laboratory coordinates (u,x,t). It is easy to see 
that this is equivalent to require that 

u r 2n -3 = M n, n = 2, 3, (65) 

with 

R(M n ) =2n-2, n = 2, 3, .... (66) 

Equation ( |65| ) and the condition ( |66] ) give us the correct terms in the expression of u T2n _ 3 but they do not determine 
their coefficients. 

Let us recall that the Maurer-Cartan equation associated with the group of 2 x 2 real unimodular matrices (Chern 
& Peng 1979), leads to the mKdV hierarchy - in variables x, t - in the form 

u t = U~ 1 R n+liX , (67) 

where 

R - -1. (68) 
R n for n > 1 may be generated from the recursion formula 

U^Rn+l^ = \( U ~ lR ^) 2x + ( URn ) x - ( 69 ) 

Let us renormalize the coefficients of highest order derivative term in each one of the equations of the hierarchy to 
+1 or —1 by defining 

A n , x = {-ir^u^R^ . (70) 
Hence, we write the higher-order mKdVI equations 

u t = A n _i )X n=l, 2, (71) 

with 

A n +i, x = -A ni3x -4Lu / uA n £ dn + u x 6 ny0 , n = 0, 1, (72) 

J — 00 

where S n fi is the Kronecker symbol. With the above statements we establish the following theorem 
THEOREM 

If ( u T3 ) = \ u T2n 3 ) with the invariance condition u T2n _ 3 — M n where 
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M 3 = -u 3i - 6u 2 U£ , 



(73) 



and 

R(M n )=2n-2, (74) 

then we have 

M n = A n - 1 ,e, n = l, 2, (75) 

where the A n are the n" 1 conserved densites of mKdVI. 
Proof: 

The commutativity of time derivatives and the invariance condition lead to the linearized mKdVI for M n 



M„, T3 + M nM + 6 (u 2 M„) ^ = . (76) 

Thus we will prove that 

( A n-i,£) r + (a„_i, ? )^ + 6(u 2 ^„-i,c) 5 = 0. (77) 

We proceed by induction. For n — 3, eq.([77|) gives mKdVI, and ([75|) clearly holds if we assume that the arbitrary 
function of integration in T3 that appears is zero. This is justified because each M n must be a polynomial in terms of 
(|62|). Assuming that it holds for n — 2, we have 

(4»-2,c) t + (A n _ 2 ,e) + 6 (w 2 J 4„_ 2 ,e) = . (78) 

We use ([72]), the mKdVI and the inductive hypothesis ( |78| ) to show - after some heavy algebra - that 

(A.-u) t + (A»-u) +6(ii 2 4- U ) =0. (79) 

This proves the theorem. 



VIII. CONCLUSION 



We have applied a multiple-time version of the reductive perturbation method to study the solitary-wave type 
solution and the kink-wave type solution of two families of modified Boussinesq model equations. 

In the first case we have eliminated the solitary-wave related secular producing terms through the use of the 
equations of an appropiated modified Korteweg-de Vries hierarchy. We have shown that the solitary-wave solution 
of the associated modified Boussinesq equation is given by a solitary-wave satisfying, in the slow variables, all the 
equations of the modified Korteweg-de Vries hierarchy. Accordingly, while the modified Korteweg-de Vries solitary- 
wave only depends on one slow variable T3, the solitary-wave solution of the modified Boussinesq equation can be 
thought of as depending on the infinite slow time variables. 

Hence, solitary-wave solutions of intermediate model equations, like the modified Boussinesq, contain complete 
information concerning all degrees of long-waves present in the system. 

In the second case, we have eliminated the kink-wave related secular producing terms through the use of another 
appropriated modified Korteweg-de Vries hierarchy. We have shown that the kink-wave solution of the associated 
modified Boussinesq equation can be built from the kink-wave solution of all equations of the corresponding modified 
Korteweg-de Vries hierarchy. 

The main difference between the solitary-wave case and the kink-wave case is situated in the formulae of tranfor- 
mation giving the slow-time variables t ,s in function of the laboratory time coordinate t. In the first case a formula 
exists and is given by the long-wave expansion of the linear dispersion relation of the initial equation, independently 
of the associated perturbation theory. 

In the second case such a formula does not exist and we know the r :S in funtion of t only "a posteriori", that is, by 
inspecting each order of the perturbative theory. 

In the solitary-wave case the sources of secularities at each order of the associated perturbative series are the linear 
terms only. If we linearize the perturbative series (u n = for n = 1, 2, ....) we obtain the equations of the associated 
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Fourier transform of the extended initial function (Sandri 1965) if, and only if, the definitions of t ,s like a function of 
t was done according to the development of the linear dispersion relation for long-waves. 

Hence, the requirement of a perturbative scheme free of solitary-wave related secularities and the existence of a 
compatible linear limit of the theory are completely equivalent. 

In the kink-wave case, there are two sources of secularity at each order: the linear terms and also some nonlinear 
ones. Hence the elimination of the linear terms, or the compatibility with the associated Fourier theory is not sufficient 
to obtain a regular perturbative series. We must carry out a second renormalization of the kink-wave frequency which 
gives us the right coefficients transforming t 2 „+i in a function of t. In this paper (Sec. VI) we realized these two 
renormalizations in only one step, but they can be realized separately. In this case the second renormalization of the 
frequency clearly appears like a reminiscence of the celebrated Stokes' hypothesis (Whitham 1974) on the frequency 
amplitude-dependence in water waves. 

Let us remember that obstacles to asymptotic integrability can be appear in the case of nonintcgrable systems 
(Kodama & Mikhailov 1995). Such obstacles was exibited in the case of the nonintcgrable RLW equation (Kraenkel et 
al 1996). For integrable systems, like the one consider here, the multiple scale method will be able to handle both, the 
solitary- wave and the A^-soliton related secularities since no obstacles to asymptotic integrability will be present. x 

Finaly in Sec. VII we have shown that the commutativity of time derivatives, together with the scale invariance 
requirement, lead to the mKdV hierarchy in t 5 , tj,..., for a field which satisfies mKdV in r 3 . 
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